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Abstract 

m ; 

We consider a spherically symmetric distribution of dust and show that it is possible, under general 
physically reasonable conditions, for an overdensity to evolve to an underdensity (and vice versa). We 
find the conditions under which this occurs and illustrate it on a class of regular Lemaitre-Tolman- 
Bondi (ltb) solutions. The existence of this phenomenon, if verified, would have the result that the 
topology of density contours, assumed fixed in standard structure formation theories, would have to 
change and that luminous matter would not trace the dark matter distribution so well. 

Of 

6 

£j : 1 Introduction 

CO 

In the inflationary universe paradigm, it is believed that the observed universe is very nearly flat. The 
density of baryons — which can be obtained from primordial nucleosynthesis theory — is however very 
small and this requires that most matter is non-baryonic. Traditional theories of structure formation 
| assert that baryonic matter fell into the high density peaks of dark matter and became luminous forming 
stars and galaxies. The stationary view, in which matter concentrations remain essentially fixed, may be 
thought of as being governed by a mapping which preserves extremal points of the density field. It may 
well be a good approximation if the initial density field is simply amplified by gravitational processing, but 
when the matter content of the pre- and post-decoupling epochs is viewed from a hydrodynamical point 
of view as a fluid in high-temperature plasma or quasi-plasma state, one would expect shock waves and 
other spatial gradients to exist (even if their amplitudes were small). Indeed, large scale inhomogeneities 
and flows have been shown to be a pervasive influence on the behaviour of the universe on scales of up 
to 100 Mpc. (Cf. for example fl3|, ||[ |[). The now undisputed existence of large-scale cosmic flows (on 
the scale of 15000 ferns -1 ) as has been reported by various authors [16], lends more credence to the idea 



that perhaps the stationary approximation, used ubiquitously in structure formation, is not as good as 
assumed. The bulk flow reading of 700 db UOkms^ 1 found for all Abell clusters with redshifts less than 
15000fcms -1 strongly excludes any of the popular models with Gaussian initial conditions. 

In this context, the Lemaitre-Tolman-Bondi (ltb) (l^, ^] universe is interesting as one may 
analytically study the evolution of spherically symmetric inhomogeneities. The discovery of large scale 
voids and walls in the eighties sparked interest in the LTB model as a means of investigating these, and 
other similar, structures (for example (l^, [p], 0]). 
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The nonlinear effects of large scale clumps (for example [[T^]) and voids ((!], [lTj) on the production 
of anisotropies in the CMB has been studied using ltb models numerically. The results have been that 
a large part of the temperature anisotropies in the background radiation (the dipole component) may 
be completely due to large scale structures, but leave open the origin of other sources (for example 
quadrupole) as truly cosmological. Also worthy of mention is the work done by Lake and Pirn (9|, |l4|]. 

These studies concentrated on the description and feasibility of spherical inhomogeneities, and 
were not too concerned with determining under what conditions structures could change radically with 
evolution. Here we intend to initiate analytical studies on this topic. 

At the centre of symmetry of an ltb universe, we must generically have a position of extreme 
density. Thus at the centre it is not feasible to study the question of density waves, per se, since a wave 
is defined by the fact that a maximum (or minimum) moves at some velocity away from the worldline. 
But at the centre we can ask the question: 1 under what conditions will a density maximum evolve into a 
density minimum or vice versa? 1 . This is a first step towards a study of cosmic flows in this model, since 
if this question can be answered in the affirmative, then it would naturally follow that in some region 
around the centre over the time elapsed a maximum (or minimum) has to be traveling away from the 
centre. 

If physical, this would raise questions about the validity of the standard model of structure forma- 
tion. It is particularly important to COBE analyses where the data (for example, hot spots and cold spots 
corresponding to under- and over-densities respectively) on the last scattering surface is 'transferred' to 
the current epoch by use of a function which does not assume that the peaks in the matter distribution 
may change to troughs. 



2 Preliminaries and Programme 

We are interested in whether the profile of a density inhomogeneity can change significantly with evolution. 
Specifically, we want to know whether a central maximum in density can evolve into a central minimum, 
or vice versa. For our investigation we use the simplest inhomogeneous cosmological solution to the 
Einstein field equations, the LTB model. 

This universe model is spherically symmetric, but in general radially inhomogeneous. Space-time 
is described by a four-dimensional continuum filled by an irrotational perfect fluid with a dust equation 
of state. We may choose the natural coordinate system labelled by {x a }^ =0 = {t, r, 9, cp} suggested by 

the spherical symmetry. The coordinates are assumed to be comoving with the particles. This allows a 

dx a 

definition of a fluid velocity u a = such that u a = 5q and u a u a = — 1, which would mean that time 

coordinate t is also cosmic time. For an ideal fluid with mass density p and vanishing pressure (dust), the 
energy-momentum tensor has the form T ab = pu a u b . The conservation of energy-momentum T ab ] f ) = 
confines the dust to geodesies and also implies that the mass of any portion of the fluid is conserved 
through the twice contracted Bianchi Identities . 

The metric in synchronous comoving coordinates can be written as 

ds 2 = -dt 2 + } R ^ dr 2 + R 2 {d8 2 + sin 2 6d<j) 2 ) (1) 
1 + 2E 

where R = R(t, r, ) acts as a transverse scale factor for individual comoving particles, and E = E(r) is 
an arbitrary function of the integration which has a dynamic as well as a metric geometric role. R is also 
the areal radius, that is AttR 2 describes the surface area of the sphere at comoving radius r at any time 
t and thus R(t, r) > 0. 

The expression for the invariant energy density, p = p(t,r), is obtained from the tt field equation: 

2M' 



where M = M(r) is another arbitrary function. 



The rr, 96 and 0</> components of the EFE reduce to the single equation of motion 



1 (R\ 2 _ M E 

2\R) ~ R* + R2 ■ W 

We define a 'scale radius', p(r), and 'scale time', q(r), for non-parabolic models as follows: 

P(r) = — (4) 

and 

q{r) = (5) 

which may be viewed as alternative variables to M and E. We (for convenience) may sometimes mix 
these four variables in the equations below. The ' + ' is applicable in hyperbolic models and the ' — ' is 
used for elliptic models. In a re-collapsing model, the areal radius at maximum expansion is given by 
p(r) and the time from creation to destruction is irq(r). 

The Friedmann-like equation can be solved (see for example @]) in terms of parameter ij = rj(t, r). 
For a non-empty universe, there are three solutions to (||): 

p 2{t-t B ) 

R = tt'Po, f = (6) 

2 q 

where 

cosh(ry) — 1, f sinh(r/) — rj, E > 

(l/2)r? 2 , i=\ (l/% 3 , E = (7) 

1 — cos(??), I r] — sm(rj), E<0 



where t B = ts{r) is a third arbitrary function. The solutions (gj7p nave the same evolution as the 
corresponding flrw dust solutions , but with spatially variable M, E and ts- In contrast to the flrw 
models, however, it is quite possible to have all three types of evolution in the same model 0. In the 
homogeneous flrw case — p = pit) only - the requirement that rj be independent of r at all times in (^, 
0) implies that 

t B = constant, M oc |£| 3/2 . (8) 

The function ts (r) is the 'bangtime function'. Individual shells of matter need not all emanate 
from one single bang event, but originate at different times as determined by ts(r). The gradient of ts 
generates the decaying modes of the perturbation to an flrw background |18| . 

M(r) is the effective gravitational mass within r. 

The local geometry is determined by E(r), as is evident from its appearance in the metric, and 
in fact this function determines the 'embedding angle' @. Also if we compare equation (||) with the 
Newtonian analogue of a dust cloud we see that E also acts as an energy potential; that is, locally 
hyperbolic, parabolic and elliptic regions occur when E{r) > 0, E(r) = and E{r) < respectively. Its 



gradient generates the growing modes of the perturbation to an flrw background [18|. 

Our method for this investigation is straightforward. We require that the density be smooth 
through the origin of our coordinate system. Thus the spatial gradient of the density has to vanish at 
R(t,r = 0) for all t. This would then impose certain restrictions on the three arbitrary functions M(r), 
E(r) and ts(r) and their derivatives for this density profile to hold. The change in concavity of the 
density profile at the origin is determined by the sign of the second radial derivative of the density at 
that point. 

For the required density profile we need expressions for the spatial gradient R' and second and 
third radial derivatives, R" and R!" respectively, explicitly as a sum of a product of functions of r and 
functions of rj. The full expressions are somewhat nasty-looking expressions and are not easily understood 
without detailed analysis, so we have merely recorded them in appendix |A[ 



2.1 Restrictions on the Arbitrary Functions 



2.1.1 Shell Crossings 

We will impose regularity conditions on the spacetime; excluding shell crossings in particular. Loosely 
stated, a shell crossing occurs when an inner spherical shell of matter moves faster than an outer shell and 
eventually bursts through. A locus of points is formed where R' = and R ^ OQ. Since the Kretschmann 
scalar K = R abcd R abcd diver ges, one may consider this to be a 'true' singularity^. In contrast to other 
studies which utilised the high-density regions created by shell crossings as generators of large-scale 
structure, we require the spacetime to be regular and thus seek to exclude shell crossings. The necessary 
and sufficient conditions under which shell crossings do not occur were derived by Hellaby and Lake ||. 



2.1.2 Behaviour at the Origin 

An origin occurs at r = when R(t, r = 0) = for all t. On any constant t surface away from the bang 
or crunch, we require that 

(a) the density p be finite, positive, and non-zero, 

M' 

- — ► Kpo(t) = const G (0, oo) (9) 



R 2 R 

(b) the Kretschmann scalar be finite 

r , 48M 2 32MM' 12(M') 2 

K = -^g- + + J^J^ -> *o(*) = const G (-oo, oo) (10) 

and (c) the evolution at r = not be different from its neighbourhood, so that (t — ts), <t>o(v) an d 

go smoothly to a finite limit in (0, oo) as r —> 0. Equation ([]) then gives us the following behaviour of 

the arbitrary functions near the origin 

R R(±E) . . _ rn . M , . 

— = * S Q {t) = const G (0, oo), q = ^ ±2E y /2 90 = const G (0, oo) (11) 

If we assume that E(r) and M(r) are analytic at r = 0, so that they can be approximated by 
polynomials in r, then we can further deduce that, as R — » 0, 

E oc ii 2 -> 0, M oc i? 3 -> (12) 

and similarly 

RocR^O (13) 
Although M'/M & both go as l/R, the foregoing gives 



q' _ fM' 3E' 
~q~ ~ \W ~ ~2E 



constant or 0, (14) 



Thus we have an FLRW-like origin] 3 ]. 



2.1.3 The Smooth Central Density Criteria 

The fractional spatial gradient of the density is obtained by differentiating the density with respect to 
R on a constant t-slice. Since R is a physically invariant quantity - the areal radius - this will give 
us results which are not coordinate dependent. We can take a slice in time in a natural way since the 
coordinate time t is also proper time for comoving dust in a synchronous metric and so also physically 

1 Regular maxima in the spatial sections also have R' = 0, but are not shell crossings fl. 

2 Other opinions are that these are non-physical in the sense that they merely indicate the impropriety of extending a 
simplified fluid description too far. 

3 If the density were allowed to approach zero at the origin, other limiting behaviours of E & M would be possible. 



invariant. Furthermore, since the 3 arbitrary functions M, E & t R all have physical interpretations, they 
are invariant too. 



Now the transformation between (t, r) and (t, R) where R = R(t, r) obeys 

dr 



Or 
dR 



r! = 1, 



Or 
dR 



R+ ,.. 
t dt 







R 



so for any function of r, say F(r) = F(r(R,t)), we define 

(IF dr _ /•" 

So the energy density on a hypersurface of constant time, equation (g), can be written as 

2d R M 



OF 
~dR 



8vrp 



R 2 



From the above equation we find the fractional spatial gradient of the density to be 

d R p d RR M 2 



where we have determined that 



d RR M 



d 2 M 



OR 2 



d R M R 



(M"-{d R M)R")/{R!f. 



(15) 



(16) 



(17) 



(18) 



We require the density to be finite and its gradient to vanish at the origin. Thus, for the required density 
profile we must have 



9rP 



P 



1 



r=0 



(R 



A2 



(10 

R 



M' 



0. 



(19) 



r=0 



We find that this gives us four conditions on the arbitrary functions and their derivatives. For the 
details in the hyperbolic case, see appendix ||. 

Requiring the density to be C 1 at the origin also implies that the bangtime function ts{r) must 
be at least C 1 at the origin. For the details see appendix |B|. 



2.2 Evolution of the Second Radial Derivative of the Density 

To answer the question raised in the introduction, we need to see what happens to the second radial 
derivative of the density. We can obtain an expression for this quantity by differentiating equation (|i"6| ) 
twice with respect to R on a surface of constant time. 

d RR p 2 (3 2 8 RR M\ | d RRR M 



p R\R d R M J d R M 
where d RR M is given by ( |i~8| ) and d RRR M is defined as 

1 :(M»'- ( a R M)iT> 



dRRR M 



dR 3 



(R') 3 " v n ' ' [R>) 2 



A more explicit form of the above for the hyperbolic case can be obtained by substituting for R', R" and 
R'" into equation (p0|). Again, this is a most unpleasant-looking expression and not easily assimilated. 
Appendix [C] contains the result we get after the assumption of a flat central density has been included. 



3 Density Profile Inversion: Existence Proof 



A change in the density contrast will depend on whether factors on the right hand side of equation 
(|56|) change sign with evolution, or terms of different sign become dominant. In order to simplify, we 
have assumed that all three arbitrary functions have polynomial behaviour at the origin. In addition, 
we imposed the smooth origin condition fll4|), which effectively says that in some neighbourhood of the 
origin, the spacetime is tangent to a homogeneous model. We now single out the dominant functions of 
rj for early and late times in hyperbolic models, which leads to the following digestible expressions. 

For early times, rj — > 0, 
9rrp 



lim 



+ 




(21) 



with R' at early times given by 



For late times we note that sinhr/ 
Therefore for large rj, and at the origin, 



R'\ 



1 M' 



rj. 



12 E 

cosh 77, cosh 77 — 1 



cosh 77 where cosh 77 — > e 71 / 2 as 77 ^ co. 



lim 

77— >oo 
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where R' at late times is given by 



16 M' E' 4 M 
¥17 ~E + 9M 1 



(22) 



R'\ 



1 MS' 



cosh rj . 



To what extent can we tailor the evolution of p by choosing the 3 arbitrary ltb functions? We 
recall that there are four restrictions on the nine quantities E, E', E", M, M', M" , tR, t' B and t" B given 
by equations (|48|)-(|5l|) ensuring a flat central density. This leaves us with the freedom to fix five of the 
above at r = 0. In addition there are also the conditions for a regular origin (section |2.1.2| ) and those for 
no shell crossings 0. However, the former conditions do not provide any additional choice restrictions 
once we have specified a flat density at the origin, and the latter are inequality constraints which only 
limit the range of choice, so are not as severe as the others. 

We see in the early time expression © there appear terms involving the three derivatives of ts 
which do not occur in the late time expression (p2j). This allows us to fix the early time behaviour of 
the density. And, likewise, the three derivatives of E occur at late times but not at early times for the 
relative change in concavity of the density. Thus, in principle, we should be able to independently fix 
the late time behaviour as well due to this freedom. In effect, we have sufficient freedom to model the 
density as being overdense initially, and subsequently evolving to an underdense state, or indeed, vice 
versa. Moreover, it is conceivable that the middle time behaviour could be separately specified, as there 
are still the derivatives of M to play with. 



4 Specific Models 

We will consider an initial overdensity changing to an underdensity. So, at early times at the origin, we 
want the concavity to be negative and at late times, positive. We consider 'exact perturbations'^ of an 

4 This is just a mathematical device. No averaging or matching procedure to define a background flrw model has been 
employed — there is no 'gauge problem' in the sense of cosmological perturbations relative to a background. 



flrw model in the following way: 



M(r) = M r 3 (l + a(r)), a(0)=0; (23) 

2E(r) = r 2 (l + /3(r)), 0(0) =0; (24) 

t B {r) = 7 (r), 7 (0) =0. (25) 

These ensure the origin conditions of section |2.1.2| are satisfied. However, they do not necessarily satisfy 
the restrictions imposed by fl4"8|)-(|5l"]). 

We are also preventing shell-crossing singularities from interfering. These conditions are, for 
t > t B , R' > and hyperbolic models ER/M > 0, 

t' B < , E' > and M' > . 

For perturbation functions of the form 

a{r) = Ar a , (26) 

(3(r) = Br b , (27) 

t B (r) = Cr c , (28) 

the requirement of no shell-crossings leads to the following restrictions on the constants A, B and C: 

t' B < C < (29) 

M' > ^ ^ > ^- — (30) 

~ (3 + a)r a y ' 

E' > ^ S > — , . (31) 

(2 + b) r b y J 



The smooth central density criteria — that is, in this case, ( |49|) and ( pT| ) for ts(r) or ( |4q ) and ( |50| ) for 
M(r) and -E(r) — impose the following restrictions on A, B and C (Tables |2| and |l|). For simplicity, we 
will only investigate models where a, b and c are (positive) natural numbers. 



Table 1 goes here 



Table 2 goes here 



Choosing a bangtime function t B = Cr 2 , substitution into equation (|2l| ) shows that the relative 
concavity of the density at early times can be fixed as negative by choosing C negative, making t B & 
decreasing function; in fact, 

160C 



,. d RR p 
nm 



r=0 Ml 



This automatically satisfies the first requirement for no shell crossings to occur as well. A choice of t B 

9rr p 

which is of higher power gives = at the origin at early times. Choosing t B as a linear function 

P 

results in a vanishing bangtime perturbation as can be seen from Table W. 

We use equation (|22] ) to determine the late time behaviour. The results after use of Table | arc 
tabulated below (Table S). 



Table 3 goes here 



Clearly there are a wide variety of models which can change concavity at the origin and which 
also have no shell-crossing singularities. 

We will illustrate the phenomenon on a model which has quadratic perturbation functions — that 
is; a, b and c are all equal to two. We choose A = 1 x 10 2 , B = 1 X 1CP 6 and C = —3 x 10 -8 . Since we 
are only interested in qualitative results, we may put Mq = 1. 

The density profile this specifies is plotted for a sequence of cosmic time (t) values in figures 1-4. 
The units (cosmological time, length, mass and density units) are converted as follows 

1 ctu = 2.005 xlO 9 yrs 

1 clu = 6.146 xlO 8 pc 

1 emu = 1.285 xlO 22 M 

1 cmu/clu 3 = 3.746 xl0~ 27 g/cc . 



Figures 1 - 4 go here 



5 Implications and Discussion 

We have shown that, for the simplest inhomogeneous cosmologies — the LTB models for spherically 
symmetric dust — a change from central density maximum to central density minimum (or vice-versa) 
during the evolution of the inhomogeneity is entirely possible, and a numerical example was presented. 
Indeed, given that the early and late time limits of the concavity of the central density profile depend on 
separate arbitrary functions that have no necessary connection, it would be surprising if profile inversions 
were not common. The models considered are completely physically reasonable for post-decoupling 
inhomogeneities . 

Perhaps the most important implications of this work derive from the existence proof of the 
possibility of density profile inversion and density waves^] in the LTB model. In the real universe, which 
is much more complex than this model, we expect waves to be generic The crucial element in our 
investigation is the importance of the bangtime function ts and its derivatives at early times. We may 
recall from section Q that t' B generates the decaying mode and E' the growing mode to RW perturbations. 
The overdensity occurs at early times because we choose ts in such a way that the result is an overdensity 
and, in a similar fashion, the underdensity occurs at a late time because we choose E such that it gives 
that particular type of density profile. Perhaps the reason why the effect obtained here has not been 
discussed before is because most studies consider linearised perturbations which have the ultimate effect 
of neglecting the decaying mode. 

The most common model of structure formation assumes Cold Dark Matter (cdm) with a Harrison- 
Zel'dovich spectrum of initial perturbations. Observations indicate that cdm predictions on large scales 
and small scales are incompatible. In particular, standard cdm has trouble reproducing the large velocity 
dispersion of luminous matter from the stationary standpoint. Realising the fairly universal failings of 
standard structure formation theories to explain bulk flow statistics, one might argue that there is some 
fundamental assumption that must be re-evaluated unless there is a radically different process responsible 
for structure in the universe. It seems natural to ask if our results might go at least some way in solving 
these problems. 

In the linear theory of structure formation, the topology of density contour surfaces does not 
change. No links are formed and no chains are broken — the genus of the surface is unchanged since 
the process is continuous. When nonlinearity is important, the genus of the surfaces evolves as clumps 
and bubbles form. Even though the statistics today may be non-Gaussian, their structure today will 
vary depending on the Gaussianity of the initial distribution. However, the way that this occurs will be 
different if the density profile inverts and if density waves are present, since no longer will overdense regions 
simply grow monotonically. There will be an interaction of spatial and temporal density gradients. This 

5 as has been numerically discovered previously in many studies (in LTB and related models) on large scale structures, 
mentioned in the introduction. 



means that density waves must be included if a correct interpretation of topological studies of structure 
formation is to be obtained. The change to the topology of the constant density contours, comparing 
density waves and no-density waves scenarios, should be examined. 

It should be clarified that the density waves discussed here are due to motion of the density 
maximum through the comoving frame, so that a galaxy that is in the density peak at one time may be 
outside it at a later or earlier time. This effect may be superimposed on the galaxy flow. 

A direct effect of this work is its implications for the transfer function used ubiquitously in standard 
structure formation theories whereby luminous matter congregates in the peaks of the underlying mass 
distributionP]. These peaks do not move; in the sense that they remain attached to the same world line as 
time evolves. The only change that happens is the infall of matter about these peaks so that the density 
contrast increases. There is spatial flow of matter, but the spatial distribution of extrema of the initial 
density field remains invariant. This invariance is broken when density profile inversions occur. The effect 
is to (amongst other things) change the form of the transfer function. We could reasonably speculate 
that the transfer function becoming more complicated may perhaps allow one to take a standard scale 
invariant spectrum and fit it to small, large and intermediate constraints. 

Another way of viewing this is that we may not be able to rely on luminous matter being an 
accurate tracer of total cosmic density, since the density peaks that triggered galaxy formation may have 
moved on. Similarly, the velocity imparted to forming galaxies may no longer be that of the unseen 
matter component. The result of the gravitational interaction of the luminous and dark components 
may be observed flows towards regions which do not seem to be density concentrations. However, such 
two-component effects are beyond the present study. 
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A Spatial Derivatives of Areal Radius 



We find 

dr 2 2 dr\ \drj J dr 

where 0o and £ are given by ([?]), and u(r) has been defined as 



After some manipulation we obtain 



„ w M' E' 
u=(M' = ---. (32) 



— = - V( j) 2 + - U( fi (33) 

ar 2 g 2 

where 0i to 0g are given in (^)-(]46|), and v(r) is defined by 

, IV 3E' M' . A ^ 

In - = — - — . 34 
qj 2E M y ' 

We proceed in a similar fashion to obtain an expression for the second radial derivative: 

-5-9- = o V ^5 - - *B «04 H o ^3 + 77 (« + 2uv)(p2 

or z 2 q q z 2 

P t B ' (w + u) 0i + J («' + u 2 )0 o (35) 



9 



where 



And for the third derivative 



W {r)=[ln*t B >) =^ + %- (36) 



d 3 R P 3 , P , / 2 , , 2p(t B ') 2 /t B 'N ; 

-5-3- = 7T V 09 - - *B ^ 08 H 5 ^0 7 - 4p 

+ (™ + v')4> 5 - ^ i B ' («' + uv + uw)0 4 

+ 5 w 03 + - (3u f + 3m> + 3u u + v )02 

q z 2 

- — t B '(v' + 2w' + 4u' + 2w 2 + 2uw + 2u 2 + uv - wv)^ 
2q 

+ |(u 3 + 3uu' W)0o- (37) 

The above derivatives of R have been expressed in terms of u, v and w because if written in terms of 
M, E and t B the expressions become a bit messy and are not very useful in that form at this stage. 
Quantities determined later will be expressed in terms of the latter variables when appropriate. 

The various functions of rj used above are 
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B The smooth central density criteria 

We want the density at the origin to be flat at all times. 

We substitute R' and R" into equation (|i~9| ) to obtain restrictions on the arbitrary functions E(r), 
M(r) and t B (r) for (H) to hold. 

We find that 




(47) 



must vanish at the origin. Here the 0j are all functions of 77 and are defined along with u(r), v(r) and 
w(r) in appendix [A]. Since the functions of parameter time 77 are linearly independent of each other, it 
follows that each of the terms in equation (E?]) must vanish separately. 



From the first term we get 



whilst the third gives 
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Note that the constraint arising from the second term is satisfied if the first (k 
are. The fourth term combined with equation (^8|) expands to 
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and third (|49T) constraints 
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The fifth term combined with the second in equation 

= 0, B 
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yields 

M" _ 
M 7 ~ 
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(50) 



(51) 



The last term produces a condition equivalent to equation (48) combined with (|50[). 

We show now that the requirement of having the density smooth through the origin, leading to 
the constraints in appendix [Bl implies that the bangtime function t_g(r) must have vanishing first spatial 
derivative at the origin. We prove this using the coordinate choice R! ~ 1 => R ~ r. As before, we can 
take the origin to be at r = without loss of generality. In these coordinates, the assumption of analytic 
arbitrary functions near r = gives p = M/(±E) ~ r 1 



q = M/(±£7)(3/2) 



r°, since M ~ r 3 and 



E ~ r . The relation which is of importance to us here is equation (|5l|) . It says that 



In this expression, we must have 

' M" M' 



Therefore, either t' B ~ r -4 to make B zero, or B = (M" /M' — 2M' /M — t" B /t' B ) diverges as 1/r or faster. 
In the former case we get ts ~ r~ 3 , which is not reasonable — either the universe is infinitely old at the 
origin only, or it will not emerge from the bang for an infinite time. In the latter case p{M" /M 1 — 2M' /M) 
is constant, so we require 

. ( M" M'\ 

tBP {~W~ 2 Jl)^ and ptB (54) 

which, near r = 0, implies 

t B ~r c , Ol. (55) 



C The Relative Concavity of the Density 



With repeated application of equations (f4§|)-(|5"l|) ensuring a smooth central density (in particular, taking 
the bangtime derivative to be vanishing at the origin); and using the variables defined by equations (|j), 
(|[), (|34]), (|50|) and (51), we find a greatly expanded form of equation (pC|) . 
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D Tables 



c = 1 


C = 


c > 2 


no restrictions 



Table 1: Restrictions imposed on the perturbation ig by the requirement of a C central density 
(VcGN + ). 





a = 1 


a > 2 


6 = 1 


A = B 


5 = 


6 > 2 


A = 


no restrictions 



Table 2: Restrictions imposed on the perturbations to E and M by the requirement of a flat central 
density (Va, 6 € 
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2,4 - 35 
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6 > 3 
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Table 3: The sign of (<9r_r p) / p at the origin at late times. 



E Figures 



Figure 1: The density profile on a worldline at an early time (t = 1 x 1CP 6 ctu ~ 2000 yrs). In these 
figures, the core is taken to be the value of R at the comoving radius r = 0.04 and they all use base 10 
logs. At this time the value corresponds to an overdensity of about 1.7 kpc in diameter which corresponds 
to the size of a small galaxy today. The units in all the figures are cosmological. The value of the density 
at the origin is po « 2.0 x 10~ 16 g/cc. 



Figure 2: The density profile on a worldline at a later time. This diagram and the next one illustrate the 
change in concavity at the centre, which occurs when the universe was ~ 2 x 10 4 years old. R ~ 7.9 kpc 
and po « 2.0 x 10~ 18 g/cc. 



Figure 3: The density profile on a worldline at a still later time which, when compared to the previous 
figure, illustrates the movement of the maximum away from the centre. This is as we expected and shows 
that a density wave must exist near the origin if the profile inverts on the central worldline. t ~ 4 x 10 4 yrs, 
R w 12.5 kpc and p « 5 x 1CT 19 g/cc. 



Figure 4: The density profile on the worldline today (t « 20 Gyrs). This corresponds to a void (albeit 
one with a rather elongated wall) with a diameter of approximately 100 Mpc. po ~ 5 x 10~ 31 g/cc and 
the maximum density on this diagram is p m f» 10~ 35 cmu/clu 3 rs 1.2 x 10 -30 g/cc. 



